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1. Introduction 

One striking property of a granular system is the ability to form a pile with a non- 
zero angle with the horizontal. This angle, called the angle of repose (pr, is influenced 
by many factors such as the grain roughness, polydispersity [1] and confinement [2]. 
Cohesive forces [31 HJ El E] also have an important influence on the angle of repose. 
A recent paper [7] has pointed out that an electric field may lead to cohesive forces 
within a granular pile. This is an interesting issue for many practical applications 
since the cohesive forces can be almost instantaneously introduced into (removed from) 
the system by switching on (off) the electric field. Understanding the effect of an 
electric field on a granular pile is also of fundamental interest. For example, it has 
been reported that ancient sediments had a larger angle of repose than present ones. A 
possible explanation [8] is the modification of Earth's electromagnetic field as well as 
the modification of electrical charging of grains due to climatic changes. The electric 
field has also an infiuence on blowing sand and dusty phenomena on Earth, Mars, the 
Moon and asteroids [9] and on the formation of geological patterns [10]. Moreover, in 
the majority of the works related to granular fiows (see [TTl [T2] and references therein) 
the infiuence of electric forces is neglected. Although this assumption is reasonable for 
large [T3l [T^ or metallic grains [15] , it probably cannot be neglected for small plastic 
or glass beads as recently demonstrated in [16]. Imposing an external well-controlled 
electric field is a good way to measure the infiuence of the electric forces on granular 
flows. 

Mersh et al. [7] list the electrically induced interactions a packing of glass beads 
can experience and compare their relative magnitudes. These interactions are: (i) 
dipole/dipole interactions due to the polarization of the grains, (ii) electrostatic force, 
(iii) interaction between two charged grains, and (iv) electroclamping force. We refer the 
reader to the work [7] (and references therein) for details. For the majority of granular 
systems, only the electroclamping force [17] is significant. We will show below that this 
is indeed the case for our system. Let us recall that when a packing of glass beads is 
confined between two electrodes a small electric current can pass through the system 
due to the presence of absorbed moisture. The density of this current increases at the 
contacts between the grains leading to a local electric field at the surface that can be 
several orders of magnitude higher that the applied electric field. Dietz and Melcher [T7] 
propose the following semi-empirical formula: F^'-' = 0.4:15TTeod'^E^;^^{EcosaY''^, where 
d is the grain diameter, E the electric field, a the angle between the electric field and 
the normal to the contact and Egat = 3 x 10^ Volts/m is the nominal electric field 
for dielectric breakdown in air (the effect of temperature, humidity, pressure. . . on this 
quantity is neglected). 

An electric field may therefore lead to the presence of cohesive forces within granular 
systems that are physically different from the adhesive forces due to the presence of 
liquid bridges [6l H]. Furthermore, the strength of these forces is set by the strength of 
the electric field, which means that it is possible to tune cohesion in a granular system 
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by modifying the strength of the electric field. Understanding the effect of an electric 
field on the stability of a granular pile will then shed light on the stability of cohesive 
granular heap. Note that cohesion can also be induced by a magnetic field [T8l[l9l|5] but 
in this case a residual magnetization might be a problem. Moreover this is possible only 
for magnetic grains whereas an electric field can induce cohesion within any granular 
system as long as a small electric current can pass through the system due to absorbed 
moisture (see below). 

Our goal is to quantify, using precise and systematic measurements of the angle of 
repose in the presence of an electric field, the influence of this field both on the static 
and on the dynamic properties of granular systems. One of the main results is that 
an electric field allows for the precise control of the angle of repose of the heap within 
a range that depends on the width of the pile. It is also possible to control the start 
or the stop of a flow on the pile by switching the electric field off or on. A theoretical 
model capturing the effect of the electrically induced cohesion as well as the effect of 
the confinement is presented. The effect of an electric field on the dynamic properties 
of the flow is found to be weak. 

2. Experimental set-up 

Experiments are performed using a quasi-two-dimensional cell. It consists of two 
vertical, fiat and parallel metallic plates (70 cm x 70 cm), as sketched in figure [T](a). 
The cell is totally closed at its bottom and on its left, so that the beads can form a 
pile. This configuration enables us to control the width W between the plates. The 
granular material we used is made of glass beads of diameter d = (500 ± 100) fim and of 
density pg = 2500 kg.m~^. Note that these grains have been used in a previous study on 
confined flow rheology [20] . In our experiments W ranges from 8d to 73d. The humidity 
(45%) and the temperature (23° C) are controlled. 

To obtain a static electric field in the cell we used a high voltage generator 
(from V^ = Volt to 3000 Volts). The corresponding electric field is, by definition, 
E = V/W. One of the metallic plates is maintained at a potential of Volts during all 
the experiments, and the potential of the other plate is freely adjustable. To form a pile 
between the two sidewalls we use the following protocol. First we impose a voltage V 
(and thus an electric field) between the sidewalls. Then, grains are continuously poured 
between the two sidewalls from a hopper. A static heap slowly forms by trapping grains 
at its top [21]. After a transient, the growth of the heap stops and the flow at its 
surface reaches a steady state. Then, the feeding of grains is cut off, and the angle of 
the remaining heap is measured once all the grains in the cell are static. In order to 
measure the angle of the heap we used two lasers as sketched in figure [1^). The first 
one is placed on a rotating axis whose height is the same as the lowest point of the cell. 
The distance between the laser and the exit of the cell, D, is fixed. The beam is always 
parallel to sidewalls and is inclined from the horizontal with an angle 6. The second 
laser is placed horizontally and its height is adjustable. The two beams are reflected on 
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Figure 1. a) Sketch of the experimental set-up. Parahel plates are made of Duralumin. 
Two lasers are used to measure the angle of repose of the heap (see text), b) Sketch 
of the pile for the calculation of the equilibrium of a thin layer. 



the pile. We adjust the height of the second laser in order to superimpose the two light 
reflections on the pile. By measuring /i and I2, respectively the height between the first 
laser and the reflection of its beam on the cell and the height of the second laser (see 
figure [I]), the angle of the heap can be calculated using equation ([T]) with a precision of 
±0.1°: 



tan I 



Dik - h 



(1) 



3. Influence of the electric field on the stability of a pile 

3.1. Electrically induced interactions 

The presence of a voltage difference on the sidewalls leads to extra forces on grains. 
Those interactions |?] are listed below and their maximal magnitude (obtained with 
the largest value of the electric field used in our experiments E = 2 x 10^ Volts/m) is 
compared to the weight of a grain. We recall that in our experiments the grains are glass 
beads (density pg = 2500 kg.m~^). Their diameter is d and their relative permittivity 

• Dipole/dipole interactions Fa due to the polarization of the grains. The ratio of Fa 
to mg is approximately 2 x 10~^ (see [7j for calculation). 

• Electrostatic force : a grain i of charge g^ immersed in the applied field can 
experience a force F^ = qiE . In order to measure the mean charge per grain Qmean, 
Ng grains are poured in a Faraday cage. In contact with one of the electrodes 
of a capacitor (whose capacity C is known) the grains charge, by influence, the 
capacitor. The measurement of the electric potential Ucap between the electrodes is 



1 X 10-1^ C, 



linked to Qmean thrOUgh : q^nean = Ucap/iNgC). We obtained Qmeai 

therefore the maximal electrostatic force Fh (for E = 2x 10^ Volts/m)) is 2 x 10~^A^. 
The ratio of F^ to mg is then 4 x 10~^. 
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Two charged grains can also interact p2j- If the charge qi and Qj of these grains 
are both equal to the saturation charge, the ratio of the corresponding force F^ to 
mg is approximately 4 x 10~^. 

The electroclamping force. Following the semi-empirical formula proposed by Dietz 
and Melcher [17]: Fj'"* = OAlbTieod"^ E^;^^{E cos a)^''^, where a is the angle between 
the electric field and the normal to the contact. If a = the ratio of F^ to mg is 
0.6. 



The only force which is not small compared to the weight of a grain is the electroclamping 
force Fd = {F^,F^,F^). This force is attractive and will increase the cohesion of the 
packing. The grain- side wall interactions can be determined by using the method of 
images. Each grain in contact with the sidewall has an image located symmetrically 
about the sidewall and with an opposite charge. Therefore the electroclamping force for 
a grain-sidewall interaction is F]'"' = 0.4:157ieod'^E°;^^E'^-'^ = k^E^-^. 



3.2. Experimental results 

For a width W between the plates varying from 8(i to 73(i, we measure the angle of 
repose without an electric field (F = 0), 0r(O) and with an electric field (pr{E) using the 
protocol described above. The results are presented in figure O The data reported in 
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Figure 2. (color online) The measured angle of repose (symbols) increases with the 
electric field. This effect, which is stronger as the gap between side walls decreases, is 
very well described by equation (|S]) (dashed lines) . 



this figure clearly demonstrate that both the confinement and the electric field have an 
effect on the angle of repose. First, without an electric field the angle of repose increases 
as W decreases. This result is in agreement with well known results |23l [21 |2l] . Second, 
in the presence of an electric field, the angle of repose can be increased by more than 
10° (for W equal to 8c? and E equal to 2.2 x 10^ Volts/m). It is therefore necessary 
to introduce cohesion and confinement effects in any model that aims to quantitatively 
predict the stability of confined and cohesive granular systems. 
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It should be pointed out that the angle of repose, for a given width W, depends on 
the current value of electric field E and not on its history. In other words, for a given 
value of E, the angle of repose will be the same if E is reached directly from Volt, 
or if several measurements for smaller value of E are made before. Starting from the 
repose angle at a given E, the angle will reach the value (/>r(0), independently of the 
previous value of E when the electric potential is switched off. By switching on the 
electric potential before the heap reaches 0^(0), the cell emptying will stop. 

3.3. Model 

We now derive a model that is able to capture the evolution of the angle of repose in 
the presence of an applied electric field. Note that this model differs significantly from 
the one presented by Robinson and Jones [25] which takes into account the effect of 
the electric field (on the grain-grain interactions as well as the grain-wall interaction) 
but not the confinement. As mentioned above, the width W between the sidewalls is 
known to have a strong influence on the angle of repose of the heap: it increases for 
decreasing W. This effect is even more important for cohesive grains |3]. Following [2] 
we can derive a simple relation between (pr and W. Let us consider a slab of material 
of width W parallel to the free surface of the pile (figure [Hb)). A surface avalanche 
starts at a given depth H when the equilibrium is broken in the direction tangent to the 
free surface of the pile. In the absence of electric field, the slab is only subjected to its 
weight as well as to the forces exerted at the two sidewalls and at the base of the slab. 
According to the Coulomb friction law, the shear forces on the two sidewalls and at the 
base of the slab are then respectively 

Fw = 2iJwNw, (2) 

and 

Fb = ^ibNb = fiBPgWLH cos 0^, (3) 

where N^ and Nb are respectively the normal force on a sidewall and the normal force 
at the base, fj,w and /i^ are respectively the effective friction coefficients at the sidewall 
and at the base, and p = upg is the density of the granular medium {u is the packing 
fraction of the medium). In the following we will assume that the effective friction 
coefficients pb and pw are constant (i.e. they do not depend on W or on y). Assuming 
that the pressure is isotropic, we have N^r = L Jq p{h)dh, where p{h) is the pressure at 
height h. The force balance along the x-direction, mg sin (pr = Fb + -FV, can be written 
as 

tan (pr = Pb + Pw TTr — TTT — ■ (4) 

WpgLH cos (pr 

From the previous equation, it can be seen that the effective friction coefficient pb is 

equal to tan^J^, the tangent of angle of the heap without walls (i.e. W — )■ oo). The 

previous equation can then be rewritten as 

tan (pr = tan (j)'^ + —^ , (5) 
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where Bm = 2/xv^/q p{h)dh/p{H). The form of this previous equation is similar to the 
one derived in [21 [26]. If we assume that the pressure is hydrostatic we have B^ = fJ-wH 
and 



tan( 



tan( 



H 

+ l^w^- 



(6) 



In figure [3] we report the angle of repose versus W without an electric field. The red 
line, corresponding to tan (pr = tan (p'i^ + ^iwH/W , is in excellent agreement with the 
experimental data. From this equation we can measure the quantity B^ = fJ'wH which 




W/d 

Figure 3. (color online) In the absence of an electric field the angle of repose depends 
strongly on the width between the sidewalls: a smaller width leads to a larger angle of 
repose. The red line corresponds to the fit tan^^ — tan0J?° + ^wH/W and the blue 
dashed line to tan 0^ = Ms = tan 0^ . 



is the characteristic length of the lateral wall effects. From our experimental results for 
E = and for the different W and using equation ([5]), we are able to determine the 
values of Bm and (f)^. We obtained Bm = 0.71Qd and 0.^ = 22.94°. The metal-glass 
static friction coefficient is between 0.5 and 0.7 pTj, therefore it is reasonable to choose 
H K. d. Moreover we will assume in the following that this height does not depend on 
the intensity of the electric field. This is reasonable as long as the cohesion force is not 
much greater than the weight of the grains. As mentioned above this is the case in our 
experiments since the maximum electric force is approximately 0.6 mg. Note that if we 
assume that equation ([6]) still holds just before the flow atop the heap stops, H can be 
seen as the flowing layer height when the flow is about to stop. In this case H ^ d 
is also reasonable for moderate cohesion. We now have to introduce the effect of the 
cohesive forces due to the electric field in the bulk and at the sidewalls. Let a^"' with 
i,j = {x, y} and a^ be respectively the stress tensors components corresponding to the 
cohesive force F^ in the bulk and the normal stress at a sidewall corresponding to the 
cohesives forces F^ exerced by the grains in contact with that sidewall. We assume that 
the stress normal to the surface is the sum of the hydrostatic pressure and the normal 
stress due to the electric force. Therefore the normal stresses at the base of the slab and 
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on a sidewall are respectively P{z = H) + cr^ and P{z) + a^ + a^ , where a^ is the 
stress tensor corresponding to the grain- wall contact. The stability condition is then 

pgHLW sin 0^ = hbWL (pgH cos 0^ + a^) +2pwL / (P{z) + (t|^ + a^j dz, 

which yields 

tan 0, = ps(l + -^^^] +f^w (§ + 2^1^^) ■ (7) 

\ pgH cos (pr J \W pgW cos (pr J 

This equation links the repose angle and the depth of the failure. Note that for a given 
(pr, two values of H can be obtained. The lowest one corresponds to a failure in the 
bulk; the highest one corresponds to a failure at the walls. This point will be addressed 
in a future paper |28] . To derive a relation between the cohesion force F^ and the stress 
tensor at the sidewalls a^ we have to determine the number of grains in contact with a 
sidewall. Recently Suzuki et al. [29] showed that the presence of a sidewall increases the 
packing fraction. In the bulk the packing fraction is close to Ubuik = p/Pg = 0.6 whereas 
it is close to Uwaii = Pwaii/Pg = 0.8 near the wall. Therefore the number of grains in 
contact with a wall of area LH is ipwaiiLH / {pgnd"^) . The corresponding normal stress 
is then 

"^ iPgndY""' 

with K = Uwall/l^bulk ~ 4/3. 

To obtain an expression for the stress components a^ in terms of electric force Fd and 
the branch vectors / joining the particle mass centers we use the formalism developed 
in [30]: (T^ = V~^ J2c&v '^iC!)F^{C)P {C) , where V is the control volume, C denotes the 
contacts in V and u{C) is the fraction of length /(C) contained in V. This requires the 
knowledge of the texture of the granular medium. For the sake of simplicity we consider 
a mean field approach where the contact network is organized as a cubic lattice. Each 
grain in the bulk has Nc = Q contacts, N-, = Ny = N^ = 2 in each direction. This 
simplification leads to a% = 26ij x 3pFf/pg7id'^, i.e. to o"|f = and to a™ = GpFy/pgiid'^. 
Equation ([7]) becomes 

tan0. = ^, + /.^ - + -— ^,^^,— ^1, (8) 



which yields 



COS -^ '^ ^"" -^ — 



^„.^,.Ata.^, = ^(2.f^ + 3F,^), (9) 

with Atan^j. = ta.n (j)r{E) — tan 0^(0). Note that in our case, Fj = F^ = k^E^''^. 

We have shown in figure [2] the comparison between the experimental data and 
the model derived above. The agreement is very good, both the effects of W and E 
are correctly taken into account. The eletroclamping force and the confinement both 
facilitate the formation of arches in the granular system. So, keeping the electric field 
E constant and increasing the gap between sidewalls should reduce the effect of the 



Electrically induced tunable cohesion in granular systems 



9 



electrical field. This is confirmed by our experimental measurements (see figure [2]). 
To test this model more carefully we also compare our data directly with equation ([9]) 
as reported in figure HI Here again the agreement is found to be very good. A small 
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Figure 4. (color online) Our model predicts that (VF/(i)Atan0r cost/)^ collapses onto 
a master curve when plotted against the electric field E. The symbols represent the 
experimental data and the dashed line represents the equation 0. 



deviation can be found for large Wjd (the model slightly underestimate the experimental 
values). This is due to the simplified structure used in our model that leads to cr|f = 0. 
For a real structure, (T|f is probably small but it is not exactly zero. 

4. Influence of the electric field on Surface flows 



Using the same apparatus, it is also possible to calculate the angle of a steady flow on 
the pile for different values of E. The hopper is now continuously filled with particles 
and its aperture precisely controls the input flow rate Q, defined as the mass of material 
entering the channel per unit of time and per unit of width [Q = mass / [time.width)]. 
When the system reaches a steady state, the input flow rate is equal to the output flow 
rate which is calculated using an electronic scale which weights the material falling out 
of the channel. The results obtained for a width W equal to 22d are presented in fig- 
ure [51 In this figure Q* is the dimensionless flow rate defined as: Q* = Q / {pgVbuikdy/gd) 
with vifu.ik = 0.6. The increase of the tangent of the angle of the free surface of the flow, 
tan(0j), with increasing Q* has already been the subject of many careful studies [3T||32]. 
Here we focus on the variation of tan(0/) with the electric field. As reported before in 
this paper, the angle of a static pile (i.e. for a zero flow rate) increases with E. The case 
of granular flows is more complex. As reported in figure [5], the influence of the electric 
field is significant only for low input flow rates. 



For large flow rates, the behavior of our system is no longer influenced by the electric 
field. The reason for this is straightforward. As mentioned before, the electroclamping 
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Figure 5. (color online) Tangent of the angle of the flow as a function of the 
dimensionless flow rate. 



force is the only non-negligible electric force in our system. This interaction exists only 
if the grains (or if the grain and the wall) are in contact. A large flow rate leads to a 
large angle for the free surface |3l] and a relatively dilute flow [33]. Thus, close to the 
free surface, the neighborhood of a given grain is dilute, the interactions between grains 
are mostly collisional and there is almost no enduring contacts. The average number of 
contacts per grain is small and the electroclamping force is thus negligible. 

5. Conclusion 

We have reported experimental results on the effect of an electric field on the angle of 
repose, (pr, of a granular heap. We show that an electric field induces cohesive forces 
between two grains in contact, and between the sidewalk and the contacting grains. 
This electrically induced cohesion increases the stability of the pile. This phenomena 
has been studied for the width between the walls ranging from 8d to 73d and can lead 
to an increase in 0^ of niore than 10°. A theoretical approach of the behavior of this 
angle is proposed, based on previous works [21 [7], including the effect of the sidewalls 
and of the electric field. The angle of a steady flow can be increased by an electric field 
but this effect is only observed for low flow rates. 
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